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ABSTRACT 

We report on the identification of a new 7 -ray-emitting narrow-line Seyfert 1 (NLSl) 
galaxy, SDSS J122222.55H-041315.7, which increases the number of known objects of this 
remarkable but rare type of active galactic nuclei (AGN) to seven. Its optical spectrum, ob¬ 
tained in the Sloan Digital Sky Survey-Baryon Oscillation Spectroscopic Survey, reveals a 
broad H j3 emission line with a width (FWHM) of 1734±104kms^*. This, along with strong 
optical Fe II multiplets [^^4570 = 0.9] and a weak [O III] A5007 emission line, makes the object 
a typical NLS1. On the other hand, the source exhibits a high radio brightness temperature, 
rapid infrared variability, and a flat X-ray spectrum extending up to ^200 keV. It is associated 
with a luminous 7 -ray source detected significantly with Fenni/hAfV. Correlated variability 
with other wavebands has not yet been tested. The spectral energy distribution can be well 
modelled by a one-zone leptonic jet model. This new member is by far the most distant 7 -ray- 
emitting NLSl, at a redshift of z = 0.966. 

Key words: galaxies: active - galaxies: jets - galaxies: nuclei - quasars: individual: SDSS 
J122222.55H-041315.7 


1 INTRODUCTION 

With the full width at half-maximum (FWHM) of the broad 
Hj 8 line less than 2000 km s^*, narrow-line Seyfert 1 galaxies 
(NLS Is*) have the strongest Fe II emission and the weakest [O III] 
emission among active galactic nuclei (AGNs). They occupy an ex¬ 
treme end of the so-called ‘Eigenvector-!’ (EVl) parameter space, 
a set of correlations between the H jS line width and other ob¬ 
servables (e.g. Boroson & Green 1992). It has been shown in pre¬ 
vious work that radio loud (RL) NLS Is are rare, only ~7% in 
fraction (e.g. Komossa et al. 2006; Zhou et al. 2006), compared to 
about 10-15 per cent in broad-line Seyfert 1 galaxies and quasars 
(e.g. Ivezic et al. 2002). The very RL NLS Is with radio loudness 
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* The term ‘NLSls’ here stands collectively for both high-luminosity 
narrow-line type I quasars and low-luminosity NLS 1 s. 


R > 100 are even rarer, accounting for only 2.5% (Komossa et al. 
2006). The reason for the rarity of RL NLSls is still unclear. 

In contrast to NLSls, the bulk of RL AGNs lies at the other 
extreme of the EVl parameter space, having broader Hj 6 and 
weaker Fe II emission (e.g. Sulentic et al. 2008). As a distinct sub¬ 
class, blazars are believed to be RL AGNs viewed with their rel¬ 
ativistic jets closely aligned to the line of sight. They are char¬ 
acterized by flat radio spectra, compact radio cores, high bright¬ 
ness temperatures and rapid variability in multi-wave bands (e.g. 
Urry & Padovani 1995). 

Since NLSls and blazars lie at the opposite extremes in 
the AGN parameter space, the discovery of their hybrids in the 
past decade (e.g., Zhou et al. 2003, 2007) has drawn consid¬ 
erable attention. By performing a comprehensive study of 23 
very RL NLSls, Yuan et al. (2008) pointed out the peculiarity 
of these objects and suggested that there exists a population of 
RL NLSls possessing relativistic jets, similar to blazars. This as¬ 
sertion was confirmed and highlighted by the detection of sev¬ 
eral objects in the GeV 7 -rays with FermilLAT (e.g., Abdo et al. 
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Figure 1. Rest-frame BOSS-SDSS spectrum of J1222+0413 after correc¬ 
tion for the Galactic extinction. The dashed line represents a single power- 
law continuum. 


2009a,b; Foschini 2011; D’Ammando et al. 2012). These discov¬ 
eries raised questions regarding the coupling of jet and accretion 
flow (Yao et al. 2015). However, further investigation of this class 
of objects in large numbers is hampered by their scarcity. There 
are only six NLSls detected in the y-ray band by FemxifLPS 
at high significance hitherto^ (PMN J0948+0022, Abdo et al. 
2009a; PKS 1502+036, IH 0323+342 and PKS 2004-447, 
Abdoetal. 2009b; SBS 0846+513, D’Ammando et al. 2012; 
FBQS J1644-I-2619, D’Ammando et al. 2015b), mostly from the 
Yuan et al. (2008) sample. 

In this Letter, we report the identification of a Fermi/LAT 
detected flat-spectrum radio quasar (FSRQ) with a new y-ray- 
emitting NLSl, SDSS J122222.55-H041315.7 (z = 0.966, hereafter 
J1222-1-0413). This object was found to be an NLSl, and the only 
NLSl associated with a FermiTLAT y-ray source, in the course of 
our ongoing program to search for RL NLSls from the Sloan Dig¬ 
ital Sky Survey-Baryon Oscillation Spectroscopic Survey (SDSS- 
BOSS; Ahn et al. 2014). J1222-I-0413 was classified as an FSRQ in 
previous surveys mainly based on its broad Mg II emission line (e.g. 
Torrealba et al. 2012). However, its NLSl nature was not known 
due to the lack of a spectrum covering the H /3 line given its rel¬ 
atively high redshift. Here we present an analysis of its optical 
spectrum obtained from SDSS-BOSS whose bandpass extending 
to ~ 10000 A covers the redshifted Hj3 region. This finding in¬ 
creases the number of the very rare y-ray-emitting NLSls to seven. 
Throughout this work a cosmology is assumed with Hq = 70 km 
s^* Mpc^*, Da = 0.73 and Dm = 0.27. 


2 OPTICAL SPECTROSCOPY 

The spectrum (Fig. 1) is first corrected for Galactic extinction 
with £(S —F)=0.016 mag (Schlafly & Finkbeiner 2011) and an 
Ri/= 3-1 extinction law, and then are transformed into the source 
rest frame at a redshift of z=0.966 provided by the SDSS spectro¬ 
scopic pipeline. The host galaxy contribution is negligible. 

We adopt a similar approach as in Dong et al. (2008) to fit 


^ Two more tentative candidates, FBQS J1102-H2239 and 
SDSS J1246-I-0238, were briefly reported by Foschini (2011). 


simultaneously the continuum, the Fe II multiplets and other emis¬ 
sion lines of the spectrum in the ranges of 2690-3100 A and 4200- 
5100 A. A single power law is used to model the sum of the thermal 
continuum emission and the contribution from non-thermal emis¬ 
sion of the jet. The optical and ultraviolet Fe II are modelled us¬ 
ing the same templates as in Dong et al. (2008) and Wang et al. 
(2009), respectively. To fit the Balmer lines, H /3 and Hy are as¬ 
sumed to have the same redshift and profile. The broad component 
of each is accounted for by either a single Lorentzian or concen¬ 
tric double-Gaussian. In addition, another Gaussian with constraint 
of FWHM<900km s^* is used for modelling the narrow compo¬ 
nent of the Balmer lines. In a first step, the [O III] XX4959, 5007 
doublet is modelled with a single Gaussian for each line. The flux 
ratio of [O III] A2,4959,5007 is fixed at the theoretical value of 1:3. 
Each of the two lines of the Mg II 2,2,2796,2803 doublet is mod¬ 
elled with two components, one Lorentzian for the broad compo¬ 
nent and one Gaussian for the narrow. The broad components of 
the Mg II doublet are set to have the same profile, with their flux 
ratio fixed at 1 . 2:1 assuming an entirely thermalized gas for sim¬ 
plicity (e.g. Laor et al. 1997). The same prescription is applied to 
the narrow components of Mg II with an additional constraint of 
FWHM<900kms-'. 


Then the models of the continuum, the Fe II and the 
emission lines are fitted simultaneously to the spectrum. We 
find that the total H /3 line is similarly well fitted with ei¬ 
ther Lorentzian+Gaussian profile or double-Gaussian+Gaussian 
profile. In the former case, the width of the broad component 
is FWHM(Hj8broad)= 1734+104 kms^* while in the latter case 
FWHM(Hj3broad)= 2264+350 km s^^. A direct measurement of 
the total Hj3 line gives FWHM(H/3total)= 1576 km s^*. These 
widths put J1222-I-0413 within the NLSl regime, near the border 
line between NLSl and normal broad-line Seyfert 1 (BLSl). In 
the following we use the Lorentzian representation of the broad 
component, because it has been shown that the broad H /3 line 
of NLSls is generally better fitted with a Lorentzian profile (e.g., 
Veron-Cetty et al. 2001; Zhou et al. 2006). In this case, the narrow 
component of Hj3 is relatively faint, with FWHM^SOOkms^'. 
The line ratio of [OIII] 2,5007 to HjStotal is ~ 0.2. The Fell is 
strong, R 4570 =Fe II A4570/Hj3total ~ 0-9 where Fe II A4570 is cal¬ 
culated by integrating from 4434 to 4684 A. These features clearly 
fulfill the conventional definition of an NLS1. 


Fig. 2 shows the residuals of the spectrum after subtract¬ 
ing the power-law continuum and the Fell model. We notice 
that an asymmetric blue wing of H /3 is present. We also no¬ 
tice that [OlIl]A5007 is relatively broad when fitted by a sin¬ 
gle Gaussian (FWHMfiil400kms^*). This may be the effect of 
an outflow in J12223-0413. It was shown in previous studies that 
NLSls often show strong blue wings and blueshifts in [Olll] (e.g. 
Komossa & Xu 2007). Thus, in a second step, we used two Gaus- 
sians to model each line of the [Olll] doublet, one for a core 
and the other for a blueshifted wing component, while the con¬ 
tinuum and other emission lines are fixed at the best-fitted results 
as above. Although the [O III] spectrum is located near the red end 
of the spectrum and has low signal-to-noise ratio (S/N) due to the 
telluric absorption, we find a broad blue wing, characterized by 
FWHM~1800kms^*, which is blueshifted by 1100kms^* from 
the narrow core component. The fitting results of major emission 
lines are summarized in Table 1. 
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Figure 2. Residuals (in grey) of the spectrum in the Mg II region (upper 
panel) and H/3+[0 III] region (lower panel) after subtracting the power-law 
continuum and the Fe II model. The best-fitting Lorentzians+Gaussians are 
displayed in black; their sum is in red. 


Table 1. Pai'ameters of the main optical emission lines. 


Line 

Model" 

FWHM* 
(km s“’) 

Flux 

(erg s“*cm“^) 

Byroad 

L 

1734 ±104 

8.1 X 10-‘5 

Byroad 

2G 

2264 ±350 

6.1 X 10-‘5 


D 

1530 

8.6 X 10^*5 

[0 in]A5007core 

G 

791 ±125 

1.0 X 10-‘5 

[0 in]A5007„ing 

G 

1841 ±671 

9.0 X 10-‘® 

Abroad 

L 

1638 ±45 

12.1 X 10-*^ 


‘^Models used for fitting the emission lines are: one Lorentzian profile (L), 
one Gaussian profile (G) or a Double-Gaussian profile (2G). ‘D’ stands for 
the direct integration over the observed line profile. 

^FWHMs are not corrected for instrumental resolution. 

3 BLAZAR-LIKE PROPERTIES 

3.1 Radio emission 

J1222+0413 was detected in several radio surveys. The spectral 
index between 1.4 and 4.85 GHz is flat, Orad = 0.3 (Sy v“. 
White & Becker 1992). The core flux density at 1.4 GHz is 0.6 Jy 
(Kharb et al. 2010), corresponding to a radio power of Pi.4ghz 
1 X 10^^WHz“^. Pushkarev & Kovalev (2012) investigated the 
very long baseline interferometry images at 2.3 and 8.6 GHz, and 
estimated a core brightness temperature of 1.43 x 10^2 and 4.34 x 
10^^ K, respectively, clearly indicating Doppler-boosted emission. 
As a blazar, J1222-1-0413 is expected to be variable in the radio 
band. We have checked the radio data in NASA/IPAC Extragalac- 
tic Database (NED) and literatures. The total flux densities varied 
in a range of ~0.66-0.8 Jy at 1.4 GHz and in a range of ~0.49- 
1.08 Jy at 5 GHz on times-cales from years to decade. The most in¬ 
tensive radio monitoring of J1222-1-0413 was performed at 15 GHz 
with cadence as short as two days, yielding the maximum variabil¬ 
ity amplitude of 24% within 405 d (206 d in the object’s rest frame, 
see Richards et al. 2011). 

To estimate the radio loudness Rsghz = 
/v(5GHz)//v(4400A), we use the 5 GHz core flux (665 mjy, 


Figure 3. WISE light curves in the wl, h’2, w3 and w4 band (from top to 
bottom panel) constructed from prohle-ht magnitudes with 1 a error. 

Laurent-Muehleisen et al. 1997) and the SDSS g-magnitude, by 
considering the k-correction with radio spectral index Orad = 0.3 
and the optical spectral index ctopt = —0.72 obtained from the 
spectral fitting (Section 2). The result is Rsghz ~ 1700. Even 
considering the radio variability, J1222-I-0413 is still a very RL 
AGN. 


3.2 Rapid infrared variability 

J1222-1-0413 was observed with the Wide-field Infrared Survey Ex¬ 
plorer (WISE) in four bands wl, >v2, w3 and w4 centred at 3.4, 4.6, 
12 and 24 flm, respectively. The variability flag in the WISE All- 
Sky Source Catalog varjiag of “9910” indicates that this source 
has a very high probability of being variable in the wl and w2 bands 
(e.g., Jiang et al. 2012). 

Its light curves are constructed from the PSF profile-fit pho¬ 
tometric magnitudes and converted to fluxes (Wright et al. 2010) 
after excluding the data whose S/N is marked as “null”. We also ex¬ 
clude those data with S/N< 10 for the wl and w2 bands, and S/N<5 
for the w3 and w4 bands, and those with reduced of the profile- 
fit photometries larger than 2. The light curves are shown in Fig. 3. 
Variability in the wl, w2 and w3 light curves are significantly de¬ 
tected within one day, with p-values of P < 0.1% using the ;|^^-test 
against the null hypothesis of no variation. This time-scale restricts 
the size of the emitting region to < 8 x 10^^ pc, much smaller than 
the scale of the torus but consistent with the size of the jet emitting 
region. The variability of w4 is not detected possibly due to the 
lower S/N and much larger photometric errors than the other three 
bands. 


3.3 High-energy emission 

Using the XMM-Newton and INTEGRAL data, de Rosa et al. 
(2008) found a hard power-law continuum with a photon index of 
r ~ 1.2 for the object’s X-ray spectrum in the broad 2-150keV 
band. A soft X-ray excess was also found below 2 keV which can 
be described by a blackbody with kT ~ 0.15 keV. 

Swift performed six observations on J 1222-1-0413 from 2007 
August to 2011 June. We reduced the data using Swift FTOOLS in 
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HEASOFT v6. 15 following the standard procedures (e.g., Yao et al. 
2015). Due to relatively poor statistics, the x-ray spectra are consis¬ 
tent with a single power law with absorption fixed at the Galactic 
value = 1.63 x 10^**cm^^, Kalberlaet al. 2005). The pho¬ 
ton indices, T =1.1-1.5, are too flat for normal radio-quiet NTS Is 
(T ~2-4, e.g., Grupeetal. 2010; Ai et al. 2011), but very simi¬ 
lar to the other y-ray-emitting NLSls (T < 2, e.g., Abdo et al. 
2009b). The 0.3-10 keV fluxes varied by a factor of 2 during the 
observations, but the correlation between the observed indices and 
the fluxes are not detected. This object is also detected with the 
Swi/r/BAT in the 14-195 keV band, as given in the Swi/r/BAT 
VOMonth Catalog (Baumgartner et al. 2013). The average hard X- 
ray spectrum is also flat, with a photon index of T =1.3. Compared 
to the previous INTEGRAL observations taken several years ago 
(de Rosa et al. 2008), the hard X-ray data show almost no signifi¬ 
cant variability (see Fig. 4). Such a flat spectrum extending to the 
hard X-ray band is characteristic of relativistic jet emission. 

In the y-ray band, J1222-1-0413 was reported to be associated 
with the y-ray source 3FGL J1222.4-I-0414 in the Third FermifLNI 
AGN Catalog (3LAC), with y-ray detection significance of 30 CT 
(Ackermann et al. 2015). The flux in the band lOOMeV-100 GeV 
is 2.9 X 10^** ergcm^^s^*, corresponding to an isotropic y-ray 
luminosity of 1.4 x 10^^ erg s^*. The photon index of the y-ray 
spectrum is 2.77 when assuming a single power-law slope, sim¬ 
ilar to those of the other y-ray-emitting NLSls which are in the 
range of 2.2-2.8 (e.g., Abdo et al. 2009b; D’Ammando et al. 2012), 
as well as those of FSRQs which are in a range of 2-3 (e.g., 
Ackermann et al. 2015). Note that y-ray sources typically are con¬ 
sidered as confirmed counterparts only, if correlated variability in 
other wavebands is detected, which is not yet seen in J1222-I-0413. 
However, we note that the hard X-ray spectrum of J1222-1-0413 is 
very flat, and its extension to the higher energies therefore is con¬ 
sistent with its y-ray detection. 

3.4 Broad Band SED 

The broad-band spectral energy distribution (SED) of J1222-1-0413, 
based on non-simultaneous data up to only 150keV (observed 
with INTEGRAL), was modelled by de Rosa et al. (2008) using 
a simple one-zone jet model. Giommi etal. (2012) used simulta¬ 
neous Planck, Swift and Fermi observations to build the SED of 
J1222-I-0413 and fitted it with a third-degree polynomial. These 
SEDs show two broad humps, characteristic of the jet emission 
from blazars. 

Here we show the SED of J1222-1-0413 compiled using non- 
simultaneous data (see Fig. 4 and its caption). We fit it with a 
one-zone leptonic jet plus a standard disc modeP using the same 
method as described in Zhang et al. (2015). We consider two cases 
for the external Compton (EC) scattering process respectively, a 
seed photon field dominated by that from the torus which is ex¬ 
pected to be present in any AGN, and by that from the broad line 
region (BLR). The energy density of the former is assumed to be 
10^^ ergcm^^ (Cleary et al. 2007) and the latter is estimated from 
the broad H j3 and Mg II line luminosities derived in Section 2 
(see Zhang et al. 2015 for details). We find that the EC/torus model 
gives a better fit to the high-energy data compared to the EC/BLR, 
and is thus considered to be the more likely model here, though the 

^ It should be noted that we adopt the standard disc to model the accretion 
flow for simplicity, which is an approximation in cases of high Eddington 
ratios. 



Figure 4. The SED of J 1222-10412 and its model. The green line represents 
the sum of all conh'ibutions from the synchrotron (syn), synchrotron-self 
Compton (SSC) and External Compton (EC) processes, and disc compo¬ 
nent. The data are taken with WISE, SDSS, Swift/GVOT, XRT, BAT and 
FermifLNI (Giommi et al. 2012). The blue trapezoid represents the X-ray 
spectrum of XMM-Newton and INTEGRAL in de Rosa et al. (2008). The 
SED of another y-ray-emitting NLSl, PMN J0948-H0022, and its models 
at different epochs taken from D’Ammando et al. (2015a, grey solid lines). 
Sun et al. (2015, grey dashed line) and Foschini et al. (2012, grey dotted 
lines) are overploted for comparison. 

latter cannot be ruled out given the non-simultaneity of the data. 
The fitted bulk Lorentz factor of the jet is Tjet = 35 for the EC/torus 
case and Fjet = 26 for the EC/BLR case. 

The best-fitting SED model (EC/torus) is shown in Eig. 4. As 
can be seen, the optical-ultraviolet band is dominated by thermal 
emission from the accretion disc, whereas the radio-to-infrared and 
the X-ray to y-ray bands are dominated by the jet synchrotron and 
EC radiation, respectively. 


4 DISCUSSION 

By analysing the BOSS spectrum of H222-I-0413, we have shown 
that the optical spectrum of this blazar is typical of an NLSl, 
very similar to PMN J0948-I-0022 - the prototype of this class 
(Zhou et al. 2003). Our finding increases the number of known 
NLSls with highly significant y-ray detection to seven. With 
2=0.966, it is the highest-redshift object of this kind detected 
so far. To estimate the black hole mass of J1222-I-0413, we 
adopt FWHM(Hj3broad)= 1734km s^* and the Hj3-based formal¬ 
ism given in Wang et al. (2009), by assuming that the BLR is 
virial as in Yuan et al. (2008). The 5100 A luminosity AL 5100 is 
estimated from the broad H j3 line luminosity using equation 5 of 
Zhou et al. (2006) to eliminate the effect of the jet contribution. We 
find Mbh ~ 1.8 X 10* (Wq (Mbh ^ 2.0 x 10* Mq if the line param¬ 
eters fitted from the double Gauissian profile are used), consistent 
with the estimation using Mg II by Sbarrato et al. (2012). Similar 
Mbh values are found when using other commonly used one-epoch 
formalisms, such as in Vestergaard & Peterson (2006). We note that 
this value is higher than the typical masses for NLSls (lO^^^M©; 
Zhou et al. 2006; Xu et al. 2012), and in fact falls between the bulk 
distributions of NLSls and blazars (e.g. Woo et al. 2005). It is also 
higher than most of the y-ray detected NLS Is, but similar to that of 
PMN J0948-I-0022 (~ 1O*M0; Zhou et al. 2003). 
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We estimate the bolometric luminosity by assuming L^oi = 
9 AL 5100 (estimated from H j3), as suggested by Kaspi et al. (2000), 
and find Lboi ~ 1.4 x 10"*^ erg s^'. The Eddington ratio of this ob¬ 
ject is A ~ 0.6, also typical of NLSls which usually have sys¬ 
tematically high A close to 1 (see Komossa 2008, and references 
therein). It is worth noting that almost all the y-ray detected NLS Is 
have black holes accreting at near the Eddington rate (Abdo et al. 
2009b). They may be analougs to the radio-bright ‘very high state’ 
in black hole X-ray binaries (e.g.. Fender et al. 2004), during which 
episodic jets are formed (see discussion in Yuan et al. 2008, and 
also Foschini et al. 2015). 

Given the close similarities in the emission-line spectra and 
the black hole mass between J1222-r0413 and PMN J0948-I-0022, 
it would be interesting to compare their SEDs. The averaged SED 
of PMN J0948-I-0022 and its model fits to the data taken at various 
epochs as given in D’Ammando et al. (2015a), Sun et al. (2015) 
and Foschini et al. (2012) are over-plotted in Fig. 4 for compari¬ 
son. As can be seen, the X-ray to y-ray emission of 11222-1-0413 
forms a single peak which can be well modelled by a strong 
EC component (this holds for both the EC/toms and EC/BLR 
cases). The SEDs of PMN 10948-1-0022 at different epochs can 
be equally well modelled by either a single or double peak, due 
to either lack of the hard X-ray detection for PMN 10948-1-0022 
and the adoption of different models of EC/torus and EC/BLR 
cases (see D’Ammando et al. 2015a and Sun et al. 2015 for de¬ 
tails). The bulk Lorentz factor of the jet for 11222-1-0413, Tjet = 35, 
according to our best-fitting model (or Tjet = 26 in the EC/BLR 
case) is similar to the value derived by D’Ammando et al. (2015a, 
Ejet = 30) in the EC/torus case for PMN 10948-1-0022, slightly 
higher than those in other y-ray detected NLSls (Tjet < 15, e.g., 
Abdo et al. 2009b; D’Ammando et al. 2012). Discoveries of more 
y-ray-emitting NLSls, especially at higher redshifts, will provide 
us with new constraints on the physics of the formation of relativis¬ 
tic jets in this intriguing class of objects and their cosmic evolution. 
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